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ABSTRACT 


Studies (Vanderburg et al. 2015; Manser et al. 2019) have shown that remnants of 
destroyed planets and debris-disk planetesimals can survive the volatile evolution of 
their host stars into white dwarfs (Villaver & Livio 2007; Duncan & Lissauer 1998), but 
detection of intact planetary bodies around white dwarfs are few (Sigurdsson et al. 2003; 
Vanderburg et al. 2020; Luhman et al. 2011; Gaensicke et al. 2019). Simulations predict 
(Madappatt et al. 2016; Mustill et al. 2018; Nordhaus & Spiegel 2013) that planets in 
Jupiter-like orbits around stars of 28M, avoid being destroyed by the strong tidal 
forces of their stellar host, but as yet there has been no observational confirmation of 
such a survivor. Here we report on the non-detection of a main-sequence lens star in 
the microlensing event MOA-2010-BLG-477Lb (Bachelet et al. 2012) using near-infrared 
observations from the Keck Observatory. We determine this system contains a 0.53 4- 0.11 
solar mass white dwarf host orbited by a 1.4 + 0.3 Jupiter mass planet with a separation 
on the plane of the sky of 2.8 + 0.5 AU, which implies a semi-major axis larger than 
this. This system is evidence that planets around white dwarfs can survive the giant 
and asymptotic giant phases of their host’s evolution, and supports the prediction that 
over half of white dwarfs are predicted to have Jovian planetary companions (Schreiber 
et al. 2019). Located at approximately 2.0kpc toward the center of our Galaxy, it likely 
represents an analog to the end stages of the Sun and Jupiter in our own Solar System. 


The microlensing event MOA-2010-BLG-477Lb was first detected by the Microlensing Observations in Astrophysics 
collaboration on 2 August 2010 (Bachelet et al. 2012). Microlensing is a technique which is sensitive to cold planets 
down to the mass of Earth (Bennett & Rhie 1996) and can probe objects around all kinds of stars, including white 
dwarfs, as unlike other detection methods it does not rely on the light coming from the host. MOA-2010-BLG-477 
is a planetary microlensing event with a planet-host mass ratio of q = (2.61 + 0.03) x 10^? and a large Einstein ring 
radius 0g = 1.26 + 0.06 mas, which implies a host star that is relatively massive or nearby. 


* Associated with Sorbonne Universités, UPMC Université Paris 6 et CNRS, UMR 7095, Institut d'Astrophysique de Paris, 98 bis bd 
Arago, 75014 Paris, France 
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Figure 1. H-band adaptive optics imaging of MOA 2010-BLG-477 from the KECK observatory a. A crop of a 
narrow-camera H-band image obtained with the NIRC2 imager in 2015 centered on MOA 2010-BLG-477 with an 8 arcsec field 
of view. b. A 0.36 arcsec zoom of the same image. The bright object in the center is the source. To the north-east (the upper 
left) is an unrelated H = 18.52 + 0.05 star 123 mas from the source. c. The field in 2018. The contours indicate the likely 
positions of a possible main sequence host (probability of 0.393, 0.865, 0.989 from light to dark blue) using constraints from 
microlensing parallax and lens-source relative proper motion. 


The target was subsequently observed with the Near Infra-Red Camera-2 (NIRC2) instrument on the Keck-II tele- 
scope on July 27, 2015, August 5, 2016 and May 23, 2018 in the H and K, near-infrared bands (1.5 — 2.3 jum) using 
laser guide star adaptive optics. Data were obtained with the wide camera (40 arcsecond field-of-view) for calibration, 
and the narrow camera (10 arcsecond field-of-view) in an attempt to resolve the source and lens. Assuming it sits 
on the main sequence, we predict a host star with mass of 0.15M < M, < 0.93M5, orbited by a planet with mass 
0.5Mjup < m, < 2.1M;jyp at a distance of 0.7 < Dr, < 2.7 kpc if all stellar types are equally likely to host a planet. 
Modelled parameters of the microlensing light curve, and the source magnitude and color (Boyajian et al. 2014), yield 
a lens-source relative proper motion of urea = 11.53 + 0.56 mas/yr. This prediction, presented here in the geocentric 
reference frame that moves with the Earth's velocity at the time of the event, allows us to estimate the future lens- 
source separation following the event's peak on 12 August 2010. 

In our Keck images we find a H = 18.52 + 0.05 star located ~123 mas to the north-east (upper-left in Fig. 1c) of 
the source. Given the large lens-source relative proper motion, a predicted separation of 57 mas means the lens star 
should be detectable under ideal observing conditions in 2015, 4.96 years after the event peak. By 2018, this separation 
will have widened to 90 mas. We find that the separation of the ~123 mas star to the north east has actually decreased 
by 11.5 mas between 2015 and 2018, instead of the expected 33 mas increase if it were the lens. The direction of this 
proper motion indicates that this object is unrelated to the lens and the planetary system. 

The source magnitude was measured as H, = 17.32 + 0.03 and K, = 17.17 + 0.04 using data from the Cerro Tololo 
Inter-American Observatory (CTIO) and the Vista Variables in the Via Lactea (VVV) survey. This brightness is com- 
patible with our NIRC2/KECK images, determined using the flux-ratio between the unrelated ~123 mas companion 
and the star at the position of the source. The sum of the two stars shown in Fig. 1c is K = 16.78+0.02 using the Keck 
data and K = 16.79 + 0.03 using the CTIO data. This is consistent with there being no excess flux within the point 
spread function (PSF) of the source star. To confirm this, we model the PSF with the addition of weak contributions 
from the wings of the unrelated companion, and find no significant structure in the residuals. Hence, with no evidence 
of a detectable lens in our Keck/NIRC2 images, the lens star must have a brightness below the detection limit. 

To determine this limit we define a detection to be above the noise at the 3o level. This corresponds to a threshold 
of H —21.1, which means that any object brighter than this value should be detectable within our Keck narrow frames. 
This limit includes an extinction correction of Ay = 0.21. The light curve of this event indicates a microlensing 
parallax signal due to the orbital motion of the Earth, but it is not fully constrained. This implies that for each light 
curve model in our Markov Chain distribution, we can determine the lens mass Mr, = (c?/AG)(0g/mg) (Gaudi 2012), 
where 7g is the microlensing parallax. We combine these light curve constraints with limits on the source and lens 
distances and velocities, together with an empirical mass-luminosity relation (Bennett et al. 2014), and find that a 
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Figure 2. H-band brightness of possible main-sequence host lenses. In white is the Einstein ring radius constraint, 
0g, derived from finite source effects in the event light curve. This constraint indicates that if the planet host star was main- 
sequence, it should be visible with Keck adaptive optics as the entire area lies above the detection threshold of our H-band 
images at H~21.1. The mass-luminosity relations for different main-sequence lens masses are shown: 1.0Mo, 0.5Mo, 0.3Mo, 
with 0.16Mo the mass lower limit derived from the mircrolens parallax, vg. Our null detection in our Keck images implies that 
the exoplanet host must be a stellar remnant, most likely a white dwarf. 


main-sequence lens star must have a brightness of H « 18.10 at 99.096 confidence. This is largely due to the light 
curve constraint on the angular Einstein ring radius, 0g, derived from the angular radius of the source star and the 
Einstein ring crossing time. 

The microlensing parallax and lens-source relative proper motion measurements constrain the location of a main 
sequence lens star to the interior of the contours in Fig. lc. The predicted brightness of a main sequence lens as a 
function of lens distance can be see in Fig. 2. Since all possible main-sequence lenses for the event are brighter than 
the Keck detection limit, and no such star is observed, the lens cannot be a main-sequence star. ‘This same analysis 
also excludes brown dwarf lenses due to an upper limit on the microlensing parallax parameter, mg < 1.03, which 
leads to an implied limit on the lens system mass of Mr, > 0.15Mo. Similarly, the lower microlensing parallax limit 
of Tg > 0.26, implies an upper mass limit of Mr, < 0.78Mo, which rules out neutron star and black hole host stars. 
With main sequence stars, brown dwarfs neutron stars, and black holes ruled out, we conclude that the lens must be 
a white dwarf. 

To estimate the properties of a white dwarf host, we use a complete sample of 130 white dwarfs within 20 pc 
(Giammichele et al. 2012), excluding unresolved double white dwarfs and double white dwarfs candidates identified 
by Toonen et al. (2017), and the Galactic model from (Bennett et al. 2014). This calculation was made under the 
assumption that all white dwarfs are equally likely to host planets. Our results are summarized in Fig. 3 and Table. 
1. We find a likely white dwarf host mass M = 0.53 + 0.11 Mo, which sits slightly below the peak of the single white 
dwarf mass distribution at 0.57-0.58 Mọ and excludes the high mass tail (Tremblay et al. 2016). This implies a Jovian 
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Table 1. System Parameters of MOA 2010-BLG-477Lb 


Parameter Units Value (+10) 20 range 


White Dwarf Lens Distance Dz (kpc) 1.99 40.35 1.31-2.69 


White Dwarf Lens Mass Mz (Mo) 0.53 +0.11 0.32-0.74 
Planet Mass Mp (Mz) 1.43 0.30 0.87-2.03 
Source Star Distance Ds (kpc) 7.81.3 5.5—10.2 
2D star-planet separation a, (AU) 2.8 3- 0.5 1.9-3.7 
3D star-planet separation a (AU) Sato 2.1-12.1 
Mass Ratio q (107?) 261+40.03 2.55-2.67 
White Dwarf Host Vz mag 26.92 + 1.80 23.38-29.44 
White Dwarf Host Jr, mag 24.93 +0.91 22.80-26.41 
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Figure 3. Physical properties of MOA-2010-BLG-477 Presented here are the (a) predicted distance to the host star/lens 
planetary system, (b) 3d star-planet separation, (c) mass of the planet, (d) mass of the host, (e) predicted host star brightness 
in V, and (f) predicted host star brightness in J. These parameters were calculated using a Bayesian analysis with priors from 
the relative proper motion of disk and bulge lenses, the lens star mass function, and velocity distributions derived from a galactic 
model as in Bennett et al. (2014) 


planet of mass mp, = 1.43 + 0.30 M; at a distance of Dr, = 1.99 + 0.35 kpc. 

With a separation from the white dwarf of a = 3.4* 3 AU (assuming a random orientation), it is likely that 
the planet MOA-2010-BLG-477Lb formed at the same time as the host star and managed to survive the post-main- 
sequence evolution. The mass loss experienced by a star on the giant and asymptotic giant branches pushes the planet 
toward a wider orbit, but tidal forces can have the opposite effect when the star expands to radii Z 1 AU (Veras 2016). 
In rare cases, the tidal effect can nearly cancel the mass loss effect, leaving a giant planet orbiting at a separation as 
small as ~ 2 AU, but this requires a fine tuning of parameters to prevent the planet from being engulfed by the star. 
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Most white dwarfs in the Galactic disk are thought to have formed from stars with initial masses of 1-2.5 Mo (Veras 
2016), and the measured separation means that to avoid tidal engulfment the mass of the progenitor star is likely to be 
< 2 Mo. Jovian planets orbiting these stars are generally thought to move to orbital separations > 5 or 6 AU around 
the remnant white dwarfs (Mustill & Villaver 2012; Nordhaus & Spiegel 2013). While this is larger than the ~ 2.8 
AU projected separation of MOA-2010-BLG-477Lb, in the region where microlensing has the highest sensitivity, the 
reduced detection probability for white dwarf planets might be compensated for by a higher intrinsic planet occurrence 
for gas giant planets around massive stars (Ghezzi et al. 2018). 

While Mercury and Venus will most likely be engulfed by the Sun when it reaches the RGB phase, Jupiter is 
expected to survive (Veras 2016). MOA-2010-BLG-477Lb is an example of the possible fate of Jupiter and supports 
the prediction that over half of white dwarfs are predicted to have Jovian planetary companions (Schreiber et al. 2019). 
It joins 29 other planets included in recent statistical sample of microlensing events (Suzuki et al. 2016). This will be 
supplemented with detections made by the Galactic Exoplanet Survey of the Nancy Grace Roman Space Telescope 
(Penny et al. 2019) aimed at determining planet occurrence rates for white dwarf hosts. 


METHODS 
OBSERVATIONS 


The microlensing event MOA 2010-BLG-477 was originally detected using the 1.8m telescope at Mt. John Obser- 
vatory in New Zealand on August 2, 2010, and subsequently observed by more than 20 telescopes (Bachelet et al. 
2012). In order to find the predicted magnitude of the source star we refer to the H-band light curve obtained by the 
uFUN 1.3 m SMARTS telescope at the Cerro Tololo Inter-American Observatory (CTIO). We calibrate the light curve 
according to data from the VVV (Vista Variables in the Via Lactea) survey (Vandorou et al. 2020). Using a single 
amplified CTIO H-band frame with an epoch corresponding to available VVV data, we cross identify between the two. 
We derive a source magnitude of HOTIO, source = 17.32 + 0.03. Given the intrinsic source color of (H — K)o = 0.07, 
and using extinction corrections derived using the OGLE extinction calculator! (Nataf et al. 2013), Ay = 0.21 and 
Ax = 0.13, we find Kerro,source = 17.17 + 0.04. 

Located at (aœ, ð) = (18^067075.47, —31°27/16.17”, J2000.0), we observed the event using the NIRC2 instrument on 
the Keck II telescope located on Mauna Kea (Hawaii) with laser guide star adaptive optics (LGSAO). Seven dithered 
H-band images were obtained with the narrow camera on 2015 July 27 (HJD — 2457230.7), 4.96 years following the 
event peak. These images have a median full width at half maximum (FWHM) of 50 mas. During this same epoch we 
obtained 14, 95 mas K,-band images with the wide camera, and 10, 50 mas H-band images with the narrow camera. 
The event was observed twice more: in 2016 (HJD — 2457605.0) with 17 H-band 50 mas narrow images and two 
H-band 90 mas wide, and in 2018 with 16 50mas K,-band narrow frames. These final 2018 observations were taken 
on 23 May 2018 (HJD = 2458262.4), 7.78 years after maximum magnification. 

The goal of these observations was to determine the composition of the source-lens blend by obtaining a refined photo- 
metric and astrometric solution. We apply flat-field and dark-field corrections using standard techniques (Blackman 
et al. 2020; Beaulieu et al. 2016). The images were stacked using SWARP (Bertin & Emmanuel 2010). We identify the 
source--lens star as the bright object in Fig. 1c. Comparing this with the OGLE-III reference image of the BLG 176.8 
field in Extended Data Fig. 1 we see that the OGLE star is a blend of four dimmer stars, the second brightest of which 
is the source--lens star. The OGLE star (number 119416) has an J-band magnitude of J = 17.446 + 0.052, which is 
consistent with the value of J = 17.443 + 0.031 from DoPhot CTIO photometry (Bachelet et al. 2012). 

We perform aperture photometry using SExtractor (Bertin & Arnouts 1996) on the wide K-band frame from 2015, 
using the narrow images to cross-calibrate the narrow frames, and find a K-band magnitude of the blend at the position 
of the source to be K Blena = 16.78 + 0.03. The magnitude is the combination of the flux from the source, the lens and 
the 123mas star located to the north-east. Fortunately, data of this field was also captured as part of the VVV survey 
while the event was still magnified. Data were obtained from the Vista telescope in JHK at the epochs: 


H, MJD — OBS = 55423.1515 
J, MJD — OBS = 55423.1576 
K, MJD — OBS = 55423.1546 


! OGLE Extinction Calculator, http: / /ogle.astrouw.edu.pl/cgi-ogle/getext.py 
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Extended Data Fig. 1. OGLE-III and Keck imaging of MOA-2010-BLG-477Lb. (a) OGLE-III image of the OGLE- 
BLG176.8 field (b) H-band image of the same field taken in 2015 with Keck/NIRC2 with the narrow camera. 


We use these data to determine whether the blend at the position of the source is consistent with the predicted source 
brightness. The image at the location of the source in the VVV image is a combination of the source, the lens, the 
123 mas star to the north-east (upper left) and the star to the north-west (upper right). The brightness of these unre- 
lated stars in the PSF are negligible when the source was amplified by 19.27 in H (where MJD — OBS = 55423.15153784 
and MJD = JD - 2400000.5). We compare the flux ratios between the star at the position of the source and the 123 mas 
companion, in the narrow Keck image and compare that to a calibrated wide Keck frame and that determined by 
CTIO/VVV. 

In our Keck wide image we crossmatch to stars in the VVV catalog (Minniti et al. 2010). Both the wide and narrow 
cameras on NIRC2 result in an image with dimensions of 1048 x 1048 pixels. The wide has a field of view (FOV) of 
40 arcseconds while the narrow has a FOV of 10 arcseconds. The narrow images were taken in sequence with a dither 
of 0.7 arcseconds. The flux ratio of between the source star and the 123 mas companion is: 


2018, K = 0.42 
2015, H = 0.33 (3 
2016, H — 0.34 


Using the flux ratio from the 2018 K data and the predicted source K brightness of K, = 17.17 + 0.04, we de- 
termine the sum of the source and the 123mas companion to be K = 16.79 + 0.04, which is compatible with 
the the Kpjeng = 16.78 + 0.03. determined from only the 2015 Keck wide co-added frame. This indicates that the 
there is no additional flux at the position of the source, and that all the photons from that object come from the source. 


LIGHT CURVE MODEL 


Since the publication of the MOA-2010-BLG-477Lb discovery paper (Bachelet et al. 2012), the MOA group has 
developed detrending methods (Bennett et al. 2012; Bond et al. 2017) that can remove systematic errors due to the 
color dependence of atmospheric refraction. We expect a significant color dependent differential refraction signal due 
to a bright main sequence star at J = 13.5 that is only 3 arc seconds from the MOA-2010-BLG-477 source star, and 
we have found a significant correction using the method of Bond (Bond et al. 2017). The best fit model parameters 
from using this new data are shown in Extended Data Table 1. Three of the parameters also apply to single lens 
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Extended Data Fig. 2. Light Curve Data and Model for microlensing event MOA-2010-BLG-477. The solid curve is the 
best fit model and the dashed grey curve is the single lens model with the same single lens parameters. The different colors 
represent different data sets from different telescopes. One sigma error bars are shown. The data sets are explained in the 
discovery paper (Bachelet et al. 2012). 


models. These are the Einstein radius crossing time, tg, the time of closest approach, to between the source and 
the lens center-of-mass, and the distance of this closest approach, uo, in units of the Einstein radius. A binary lens 
requires three additional parameters, the lens separation, s, in units of the Einstein radius, the angle between the lens 
separation and source trajectory, a, and the planet-star mass ratio, g. Events like MOA-2010-BLG-477 with caustic 
crossings require the source radius crossing time, tą, to account for finite source effects. Finally, it is also important to 
include parameters to describe the orbital motion of the lens system and the Earth-based observers. The orbital period 
of the planet is many years, so the only detectable effect of orbital motion is the relative lens velocity, $;, and 5,, in 
the directions parallel and perpendicular to the lens separation. The microlensing parallax is described by the North 
and East components of the microlensing parallax vector, mg. The inclusion of microlensing parallax introduces an 
approximate degeneracy that corresponds to a flip in the orientation of the lens system with respect to the orientation 
of the Earth's orbit. We label the solutions with the sign of the ug parameter. The uj > 0 model is favored by 
Ay? = 5.08, and it is displayed in Extended Data Fig. 5, but both models are compatible with our overall conclusions. 

One major difference between the analysis in the discovery paper (Bachelet et al. 2012) and our new analysis 
is that the Einstein ring crossing time, £g, has dropped by about 16%. This causes a similar increase in the mass 
ratio, q, because the timescale of the planetary signal has not changed. The other significant change in the model 
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Extended Data Table. 1. Best Fit Model Parameters 


This Analysis 


Discovery Paper 


parameter uo > 0 uo « 0 uo > 0 uo « 0 
tg (days) 39.523(425) 39.796(441) 46.868 47.560 
to (HJD’) . 5420.9381(3) 5420.9376(3) 5420.9382  5420.9374 
uo 0.003966(51)  —0.003937(54) 0.003430 —0.003315 
5 1.12203(52) 1.12278(51) 1.12407 1.12429 
a (rad) 0.54126(176) —0.54137(174) 0.54004 —0.53836 
q x 10° 2.590(31) 2.587(31) 2.210 2.168 
tx (days) 0.02736(18) 0.02739(18) 0.02719 0.02719 
$z 0.00039+6:60040 0-00030T3 6002 0.00223 0.00016 
Sy —0.001927 9 00746 0.002744500394 —0.00587 0.00196 
TE,N 04450 0205. — 530090 5207-— 1005432 0.0974 
TE,E =0.3165 0225 0.14670 oe  —0.1026 ^ —0.0134 
x? /dof 6420.91/6490 | 6425.97/6490 . . 


NOTE-—AII errors are given as +lo. 


parameters can be found in the microlensing parallax values. As shown in Extended Data Fig. 5a, only one component 
of the microlensing parallax vector, mg, is tightly constrained by the data. It is quite common for the component 
parallel to the direction of the Earth's acceleration at the time of the event to be measured much more precisely 
than the perpendicular component, and in this case the uncertainty in the perpendicular component is increased by 
the degeneracy between microlensing parallax and lens orbital motion. This parallel component is roughly in the the 
East-West direction for events in the Galactic bulge. However, the full mpg can be determined (Bhattacharya et al. 
2018; Bennett et al. 2020) with a measurement of the lens-source relative proper motion, Lac, although this must 
be be done in the Geocentric frame that has been used for the light curve modeling. 

While the changes in the model due to the improved MOA data are quite noticeable, they have very little effect on 
our conclusions. We have performed the full analysis using the light curve models from the discovery paper (Bachelet 
et al. 2012), and these implied very similar results. The only difference is that the old analysis shifted the tg, g toward 
more positive values by 0.1—0.2. As a result, the old light curve analysis allowed higher mass host stars. 

We note that the discovery paper suggested that it might be sensible to search for three body lens systems that could 
explain this light curve. However, unlike the case of published triple lens events (Gaudi et al. 2008; Bennett et al. 
2010; Gould et al. 2014; Bennett et al. 2016), this light curve shows no evidence of an unmodeled light curve feature 
that can be explained by an additional lens. Thus, we can expect that the only effect of an additional lens would be 
to increase the uncertainties in the properties of the lens masses in the binary lens model (Zhu et al. 2014), usually 
by very small amounts. Since our primary conclusions are based on the properties of the host star, an additional lens 
that is consistent with the data would not affect our conclusions. 


BAYESIAN ANALYSIS LIGHT CURVE MODELS WITH A GALACTIC PRIOR 


Our analysis makes use of the light curve models from the Markov Chain Monte Carlo calculations from our reanalysis 
of the MOA data with our detrending method (Bond et al. 2017), as well as other data sets from the discovery 
paper (Bachelet et al. 2012). These models incorporate the effects of the orbital motion of the Earth which are 
responsible for the microlensing parallax, as well as the orbital motion of the planet. The microlensing parallax is 
a two-dimensional vector, mg, parallel to the lens-source relative proper motion, which means the distribution of 
microlensing parallax vectors from the light curve models constrains the direction of the lens-source relative proper 
motion. The amplitude of the lens-source relative proper motion vector is determined by finite source effects in the 
microlensing light curve. However, the direction and length of this lens-source relative proper motion vectors are 
not determined in the Heliocentric reference frame that is appropriate for high angular resolution Keck follow-up 
observations. The light curve models employ a Geocentric reference frame that moves with the velocity of the Earth at 
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Extended Data Fig. 3. White Dwarf Mass-Luminosity distribution derived a sample of 130 white dwarfs from a 
homogeneous and complete sample of white dwarfs within 20pc of the Sun Giammichele et al. (2012). Two unresolved double- 
white dwarfs (DWD), eight unresolved DWD candidates and one unresolved binary white dwarf with a main-sequence companion 
have been removed from this sample (Toonen et al. 2017). 14 stars with distances > 20 pc have also been removed. The white 
dots indicate the masses and V band magnitudes of the white dwarfs in this sample, and the color distribution indicates the 
smooth Gaussian multivariate kernel-density distribution that we have used in our analysis. 


the time of closest lens-source alignment. For MOA-2010-BLG-477, this velocity is vg v.g = (—2.7933,19.5634) km s^! , 
and the transformation to the relative proper motion in the Heliocentric reference frame, rei, is given by 


NUS Trel ; (3) 
(Dong et al. 2007; Bhattacharya et al. 2018), where the relative parallax is given by Tye) = AU/Dr, — AU/Dsg. This 
transformation from Hrel,G to Mrel,H cannot be computed with light curve parameters only as mre; depends on the source 
and lens distances. Therefore, we must invoke a Galactic model to properly sample the source distance, Dg, values. 
We have used the Galactic model from Bennett et al. (2014) in our analysis. Our results are largely independent of 
the choice of Galactic model for this Bayesian analysis. The only poorly understood prior is the host mass dependence 
of the planet hosting probability, since we are considering a lens that is known to host a planet. For simplicity, we 
assume that all white dwarfs are equally likely to host a planet. Once the source distance is selected, the lens distance 
can be determined from light curve parameters, 


Mrel,H = Hrel,G T 


AU 
D, = — ~ . 4 
5 nghe + AU / Dg ( ) 
While 7g is only partially constrained by the light curve data, it is precisely determined for each model in the Markov 
Chain, so equation 4 can be used for each Markov Chain model. The Galactic model also provides weights for the 
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Comparison of best and best static models 
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Extended Data Fig. 4. Cumulative Ay? comparing a static light curve model with that including parallax and 
orbital motion. The bulk of the signal comes following the light curve peak. The parallax plus orbital motion comes primarily 
from the MOA data (Ax? = 45) and the SAAO data (Ax? = 9.0). 


density of stars at the distance to the lens, Dr, and by the probability of a star with a mass equal to the lens mass. 
This implicitly includes the assumption that all stars are equally likely to host the planet of the measured mass ratio 
and projected separation. 

Extended Data Fig. 5 shows how the Galactic model affects the distributions of the mg and Hrem vectors. Panels 
5a and 5b show the mg distribution from the light curve models in (Bachelet et al. 2012) and the implied prey 
distribution, with the help of the Dg distribution from the (Bennett et al. 2014) Galactic model. The component of 
the mg parallel to the direction of the Earth's acceleration during the event is tightly constrained. This component 
is close to the East direction. Conversely the perpendicular component, which is largely in the North-South direction 
is very weakly constrained. Panels c and d of Extended Data Fig. 5 show the result when we apply the complete 
Galactic model including the mass function for main sequence stars. This removes the light curve models with small 
Tp values, and therefore, large masses. It was these low-7p, high-mass light curve models that allowed the me vector 
to point in any direction. This was responsible for the ring distribution of Hrem, but with the high mass lens systems 
excluded, this ring is broken into two arcs to the North and South, with low-mass lens systems only allowed in the 
Northern arc. Our Keck observations have ruled out any main sequence stellar lenses in these arcs. Our model does 
assume that stars with masses above 1.1 M have left the main sequence, but such stars would also be brighter than 
the source star, which are clearly ruled out over the full ring in Extended Data Fig. 5b ring. 

With main sequence hosts ruled out, we can now turn our attention to white dwarf host stars for the MOA-2010-BLG- 
ATTLb exoplanet. We can repeat the same calculation that we did for main sequence sources. This requires a mass- 
luminosity distribution for white dwarfs. We construct such a distribution from the 20 pc sample as in (Giammichele 
et al. 2012), excluding unresolved binary white dwarfs because these are likely to have unreliable parameters. The 
resulting white dwarf mass-luminosity relation, constructed using a multivariate Gaussian kernel density-estimation is 
shown in Extended Data Fig. 3. The results of a repeat of the Bayesian analysis with the main sequence mass function 
replaced by our white dwarf mass function, from Extended Data Fig. 3, is shown in Extended Data Fig. 5 panels e and 
f. The results are quite similar to the main sequence case, but the exclusion of low-7zg, high-mass light curve models is 
now somewhat stronger because the white dwarf mass function is strongly peaked around Mz ~ 0.57 — 0.58M. This 
Bayesian analysis is also used to produce the lens system properties presented in Table 1. The contours from Extended 
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Extended Data Fig. 5. Predictions of the microlens parallax vector zz and the corresponding predicted relative 
lens-source proper motion Hre; for a main sequence and white dwarf lens. Based on a Markov-Chain Monte-Carlo (MCMC) 
analysis using Galactic model priors as in Bennett et al. (2014), the upper panels (a) and (b) show the unweighted predicted 
components of (ngen, TEE) and (Urel, HN, Hre, HE). The middle panels (c) and (d) show the weighted predictions for a main- 
sequence lens. The lower panels (e) and (f) show the weighted predictions for a white-dwarf lens. The three shades of blue from 
dark to light denote probabilities of of 0.393, 0.865, 0.989. When integrating over all parameters the limit of the 0.393 contour 
corresponds to the 1c distribution of any chosen parameter. 
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Extended Data Fig. 6. H-band adaptive optics imaging from the KECK observatory, with contours showing the 
predicted position of a white dwarf lens (analogous to Fig. 1) (a) A crop of a narrow-camera H-band image obtained 
with the NIRC2 imager in 2015 centered on MOA 2010-BLG-477 with an 8 arcsec field of view. (b) A 0.36 arcsec zoom of the 
same image.The bright object in the center is the source. To the north-east (the upper left) is an unrelated H = 18.52 + 0.05 
star 123mas from the source (star 123NE. (c) The field in 2018. The contours indicate the likely positions of the white dwarf 
host (probability of 0.393, 0.865, 0.989 from light to dark blue) using constraints from microlensing parallax and lens-source 
relative proper motion. 


Data Fig. 5d are reproduced in Extended Data Fig. 6c, where they replace the contours for the main sequence stars 
shown in Fig. 1c. 


DETECTION LIMITS AND POINT SPREAD FUNCTION DECONVOLUTION 


We determine the detection threshold of our Keck images by estimating the flux of a point source and evaluating the 
normalized cross-product with the point-spread function (PSF). We obtain this quantity for all points in the subtracted 
image and construct a map of the amplitudes. In order to obtain an estimate of the fluctuations due to the noise we 
calculate the standard deviation in this map. We consider that these fluctuations are significant at the 30 level and 
that this 3c level is our detection limit. We then convert this 3c limit to magnitude and obtain a minimum limit of 
detection of H ~ 21.1. 

While the relative proper motion indicates that the lens is distinct from the PSF of the source star, we then perform 
an analysis to determine if the star at the position of the source is consistent with a single object or if there is evidence 
of a two component system. To do this we make a numerical estimate of the PSF by stacking the brightest stars 
in the neighborhood of our target. The accurate position of each PSF is estimated by iterative Gaussian weighted 
centering. The PSFs are then interpolated, re-centered and stacked, and median solution obtained. In our H-band 
images the star of interest is quite close to the star to the north-east, which means it may receive some contribution 
from the PSF wings of its neighbour. To subtract any neighbour contributions we reconstruct the wings of the PSF of 
the more distant star as a single function of distance and subtract this weak contribution from the image. Then the 
best solution obtained is subtracted, leaving only the stars of interest. The result of this subtraction and the residuals 
are presented in Extended Data Fig. 7. 

When reconstructing a single star a common problem is to determine whether this star is a single PSF or a very 
close system of multiple PSFs. If the single PSF fit results in significant residuals, for example, it is clear that the fit of 
a binary system should be attempted. In this case constraints on the PSF can be used to find upper limits on allowed 
separations that are compatible with the data and noise. Our method to achieve this and to overcome the degeneracy 
created when the PSF components are very close are described below. 

The image as shown in Extended Data Fig. 7, Ic is the result of the convolution of the object data, Io with the PSF, 
@. Let's consider the case where the un-convolved data Jp is very narrow with respect to the PSF, which is typical in 
cases of two close-together stars. In this instance, we can write: 


ioe J d(z — uy — v)Io(u, v)dudv (5) 
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Extended Data Fig. 7. Keck Point Spread Function (PSF) fit and residuals. (a) Keck/NIRC2 H-band image from 
2018. (b) Residuals after fitting the PSF using multiple stars in the neighbourhood of the source. Both the object to the 
north-east (upper left in panel A) and the object at source position (lower right) are subtracted using this PSF fit. There is no 
structure or indication of a double star in either of the two objects. (c) The residuals from panel B normalized to the Poisson 
noise. (d) Panel A but subtracting the fitted PSF from the unrelated companion. 


Considering that the variations are small with respect to the variations of ¢, we can write $ as a local expansion in 
the local variables (u, v): 


` 06 p Pe GH $ 
Io(x,y) = moó(z, y) ML 2 By Tt ar j 4 aay "is gy? B 


where: 


mo — f Io(u, v)dudv 

mı = f Io(u, v)u dudv 
m» = f Io(u, v)v dudv 
ma = 4 f Io(u, v)u?dudv 
ma = f Io(u, v)uv dudv 


ms = 5 f Io(u, v)v?dudv 


The mı and mo coefficients represent the degrees of freedom related to the centered of the function Jp. To eliminate 
these degrees of freedom we make the center of flux of Jj coincide with the origin of the coordinate system. In this 
case, m; = 0 and mz = 0. As such we are left with an expansion with 4 basis functions, the PSF and its 3 second 
order derivatives. The moments of the function Io are simply the 3 coefficients (ma, m4, ms) normalized by the total 
flux (Note that provided that the PSF is normalized the total flux and the coefficient mo should be very similar). 
We use the numerical model of the PSF to reconstruct the derivatives up to the second order. The derivatives are 
obtained by shifting the PSF model and taking the difference with the original PSF. The value of the shift is small 
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with respect to the size of the PSF grid. We choose a value of 0.01, but also tested 0.1 and 0.001 and no significant 
changes were observed. In creating the PSF model we take an area around the object large enough to include the PSF 
wings but small enough to avoid including another object. 
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DATA AVAILABILITY 


The KECK Observatory data used in this study are freely available on the Keck Observa- 
tory Archive (https:/ /koa.ipac.caltech.edu/cgi-bin/KOA /nph-KOA login). Data from the VISTA  Vari- 
ables in the Via Lactea (VVV) survey are available on the European Southern Observatory archive 
(http: / /archive.eso.org/wdb/wdb/adp/phase3 main/form?phase3 collection- VVV &release tag—6). Data used to 
model the light curve are available from the authors on request. 


CODE AVAILABILITY 


The KECK pipeline is available on Github (https:/ /github.com/blackmanjw /KeckPipeline). The Bayesian analysis 
code of D.P. Bennett uses routines from (Press et al. 1992) which are subject to restricted availability. 
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